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SUMMARY 


This paper extends the author's prior development of a general abstract 
representation for the design sensitivities of Green's functional for linear 
structural systems to the case where the structural stiffness vanishes at an 
internal location. This situation often occurs in the optimal design of 
structures. Most optimality criteria require that optimally designed beams be 
statically determinate. For clamped-pinned beams, for example, this is 
possible only if the flexural stiffness vanishes at some intermediate location. 
The Green's function for such structures depends upon the stiffness and the 
location where it vanishes. A precise representation for Green's function's 
sensitivity to the location of vanishing stiffness is presented for beams and 
ax i symmetric plates. 


INTRODUCTION 


This paper is concerned exclusively with the linear self-adjoint differ- 
ential equation, represented in abstract form by 

Lu * T* ET u = f in fi (la) 

£ 

Here T and T are operators which are L 2 ( ) adjoints of each other, E is a 
stiffness operator which is symmetric with respect to the L 2 (ft) inner product, 
u is the response function and f is a specified disturbance. The open region 
\i e R n is bounded by 3Q. 

Appropriate mixed inhomogeneous boundary conditions are appended to 
equation (la). These are 


BY u - g on 3f^ 
B*Y*ET u = h on 3L 2 


(lb) 


where oQ-j U3f*2 = and olii H 3P 2 = <P« The operators Y and Y* map functions in 
the domain of L into functions defined on 3fii and 3Q 2 , respectively. And the 
operators B and B* map functions defined on dft-| and 3ft 2 , respectively, into 
functions defined on 3fli and 3fi 2 . Examples of the operators appearing in 
equations (la) and (lb) can be found in reference 1 for a number of specific 
applications. 


*This research was supported by the Howard University Large Space Structure 
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The stiffness operator E frequently depends upon one or more design parame- 
ters, which are collectively denoted by S. The operators T and T* are gener- 
ally differential operators which are independent of the design. The boundary 
operators may or may not be design dependent. There are two important classes 
of problems for which the boundary operators depend upon S. One class of such 
problems is usually referred to as shape optimization problems. Here, the 
boundary is the design variable, and consequently the boundary operators are 
necessarily design dependent. The other class of such problems occurs in 
structural optimization theory whenever optimality requires that the stiffness 
vanish somewhere in the interior of ft. In this case, equation (lb) must also 
include an internal boundary where certain jump and/or continuity conditions 
are specified. This latter class of problems is the primary concern of this 
paper. 


GREEN’S FUNCTION AND FUNCTIONAL 

Oden and Reddy (ref. 2) have shown the operator P, which consists of the 
spatial operator of equation (la) and the boundary operators of equation (lb), 
will be self-adjoint if the following integration by parts formula is satisfied 

(Tu.ETvL " (u ,T*ETv ) - (Yu,B*Y*ETv) - (BYu,Y*ETv) (2) 

ft ft 3ft 2 8ft x 

for every u and v in the domain of P. In equation (2), (.,.) denotes the usual 
L 2 inner product and the appended subscript the domain of integration. Thus, 
for example, (*,*K^ denotes the L 2 Oft ) inner product. In the remainder of 

this paper, it will be assumed that the operators specified in equations (la) 
and (lb) do indeed satisfy equation (2). 


The solution to equations (la) and (lb) can now 
Green 1 s function G, corresponding to the operator P, 


be obtained in terms of 

i .e. 


u - <r,G) a • (g.Y ETC) 3ni . (h,YG) 3a2 
Equation (3) may be routinely derived by noting that G(x,y) satisfies 


(3) 


T* ETG(* ,y) = 6y in ft 

and boundary conditions 

BYG( • ,y) - 0 on 3ftp 

B*Y*ETG(-,y) = 0 on jft 2 

where 6 y represents the Dirac distribution with a singularity at the location y. 
Upon taking the L 2 (ft) inner product of both sides of equation ( t»a ) with u and 
integrating the result twice by parts according to equation (2), equation (3) 
immediately follows. Several illustrations of equation (3) have been derived 
by Roach (ref. 3) for specific operator equations. 


(4a) 


(4b) 


Green's function G(x,y) is defined on the Cartesian product space ft x ft 
and is generally singular when x=y. If any of the operators appearing in 
equation (4) depend upon the design variable(s) S, then G is a functional of 
the design S. Reiss (ref. 4) recently presented a compact formula for the 
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design derivative of G when E is the only operator appearing in equations (la) 
and (lb) that depends upon the design S. For, in this case, let S and S + AS 
denote two designs and define AG by 

AG(x ,y ;S, AS) = G(x,y;S+AS) - G(x,y,S) (5) 
It immediately follows from equations (iJa), (^b) and (5) that 


T*ET AG - 

F 

in 

ft 


BY AG = 

0 

on 

9ft i 

(6) 

B*Y*ETAG = 

H 

on 

dft 2 



where 

F*- T*AET(G+AG) 

* * (7) 

H a - B* i AET(G+AG) 

A cursory comparison of equation (6) with equation (1) shows that the so- 
lution for AG is immediately specified by equations (3) and (7); thus 

AG - - (T*AET(G+AG) , G)^ 

- (bVaet(g+ag), YG) (8) 

O J62 

After applying the integration by parts formula (2), the variation AG simpli- 
fies to 


AG “ - (TG, AET(G+AG)) n (9) 

Equation (9) is an integral equation for AG. Considerable simplification re- 
sults if E is Gateaux differentiable with respect to the design. In this case, 
by restricting the design variation AS to be infinitesimal, AE is also infini- 
tesimal and equation (9) may be linearized, i.e., 

6G » - (TG, 6ETG)^ (10) 

In equation (10), the symbol A has been replaced by 6 in order to signify 
linearization. Equation (10) represents the design sensitivity of Green’s 
functional . 


SINGULAR DESIGNS 


Beams 

As stated at the outset, the primary focus of this paper is on designs 
whose stiffness vanishes somewhere in the interior of 0. For beams whose 
boundary conditions are specified by (lb), the stiffness vanishes, at most, at 
two internal locations. Let x D denote the typical location for which S(x 0 ) = 0. 
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In terras of conventional notation for beams, the internal boundary con- 
dition at x 0 is the prescription of zero moment, while the matching condition 
is zero jump in both the shear force and the response. Thus 


S ( x ) Gxx( x *y* x c>) 1 x=x ” 0 

0 

(lla) 

[[(S(x)G xx (x,y;x c ) ) x ]] x ^ Xq - 0 

(lib) 

o 

ti 

X 

!l 

X 

i — i 
i — \ 

c 

X 

>> 

X 

0 

1 1 

l — 1 

(11c) 


o 


where the subscripts denote partial differentiation with respect to the 
indicated arguments, and [[• ]] denotes the jump in the quantity within the 
double brackets. In addition, at the extremities of the beam, G must meet the 
static and kinematic boundary conditions specified by equation (4b). 

If the neighboring design S + 65 also vanishes at x D , then 6S is specifieG 
by equation (10). If, however, S + 6S vanishes at x D + 6x 0 , then 6G depends 
explicitly upon 6x 0 as well as 6S. Since the sensitivity of G with respect to 
6S is determined by equation (10), it remains to investigate the sensitivity of 
G to variations in x D . 

With x G treated as the design variable, the counterpart, to equation (5) 
becomes 


AG(x ,y ; x D , Ax c ) a G(x ,y ,x 0 +Ax 0 ) - G(x,y;x 0 ) 

which, upon linearization, simplifies to 

6G -- G x (x,y;x 0 ) 6x 0 (12) 

o 

It is important to note that G will generally have a slope discontinuity at x 0 , 
but G + 6G will have a slope discontinuity at x D + 6x 0 . It follows from 
equation (6) that 6G satisfies 

(S6G XX ) XX “0 0 < x < x G , x D < x < L (13) 


plus appropriate homogeneous boundary conditions at x=0 and L. Due to the 
shift in the internal boundary x D , care must be taken in determining the in- 
ternal matching conditions for <5G. While G satisfies equations (11a), (lib) 
and (11c), G + 5G must satisfy 


S(G+SG) XX | X -x + gx “ 0 
o o 

[ [ ( S( G+ 6G) xx ) x ] ] x „ x + (S X = 0 

o o 

t[G+6G]] X ax + 6x ~ 0 


Next, SG XX is expanded in a Taylor series about x D to get 

SG I , , = SG I + (SG ) I 6x 

XX lx + 6x xx ! x=x xx'x'x o 

O O O o 

which, by virtue of equations (11a) and (14a), becomes 


S6G 


(SG ) 

— w ' - 


6x 


XX x-x Q = ' xx'x'x=x 0 o 


(14a) 

(14b) 

(14c) 


(15a) 


322 



Similarly, by expanding (SG XX ) X and G about x « x D , and making use of equations 
(lib), (11c), (13), (l^b) and (l^c), the following jump conditions are ob- 
tained: 


[[S6G XX ) X ]] X=X = 0 (15b) 

o 

[[6G]] x=Xo = - [[G x ] ] x=Xq 6x 0 (15c) 

The sensitivity oG is completely specified by equations (13), (15a), 

(15b), (15c) and the boundary conditions at x=0 and x=L. After multiplying 
both sides of equation (13) by G and integrating the result over the domains 
0 < x < x D and x 0 < x < L, it is found that 


6G(z,y) — L[G(x,z)(S(x)6G xx (x,y) ) x x u X 

- [[G x (x,z)S(x)6G xx (x,y)]] x „ Xo 
+ [[S(x)G xx (x,z)6G x (x,y)]] x=x 

- [[(S(x)G xx (x,z) ) x 6G(x,y) ]] x= 


x=x 


(16) 


Equation (16) can be considerably simplified by making use of the jump con- 
ditions (11a, b,c) and (15a, b,c). The first term on the right hand side of 
equation (16) vanishes by virtue of equations (11c) and (15b); the third term 
also vanishes as a consequence of equation (11a). Now, substitution of 
equation (15a) into the second term, and equations (11b) and (15c) into the 
fourth term yield 

SG(y ,z) — — { E C G x ( x , z ) j j Xi . x Q(x 0 ,y ) 

o 

+ [[G x (x,y)j] x=x Q(x c ,z) } 6x 0 (17) 


where Q (x Q ,y) is the shear force at x 0 due to a unit load at y. Thus 


0(x o ,y) - “ (S(x) G xx ( x , y ) ) x | x=x q (IS) 

The design derivative of Green* s function, obtained from equations (12) and 
(18), becomes 


3G(x ,y ; x 0 )/ 3x c - “ [ LG X (x ,z ) J ^(Xq,^) 

" [ E G x ( x ,y ) ]] x - x ^ Q ( x q , z ) 


(19) 


Axisymmetric Circular Plates 

Thin isotropic elastic plates, like the elastic beams considered above, 
obey the fundamental equations (4a) and (4b), Consequently, the sensitivity of 
Green* s function with respect to changes in the plate stiffness (thickness) 
must satisfy equation (10), For simplicity, only ciruclar plates subject to 
axisymmetric loads and boundary conditions are considered. The plates may be 
full or annular, and the inner and outer boundaries of the plate will be de- 
noted by a and b, respectively. For full plates, a=0. If the stiffness of the 
plate vanishes over a circle of radius r 0 and this radius is also a design pa- 
rameter, then the sensitivity of G with respect to r 0 also must be determined. 
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At a circle of vanishing stiffness, the radial bending moment vanishes, 
and both the shear force and the response are continuous. Thus the counter- 
parts to equations (11a), (lib) and (11c) are, respectively, 

S( r ) |rG rr (r , £ ; r D ) + vGp(p,£;PQ)}|p = p — 0 

o 

[[r(S(r)(rG rr (r,c;r 0 ) + vG r (r , x, ; r 0 ) ) ) r (20) 

— S(r)(G r (r,c;r 0 ) + vrGpp(r,^;rQ))Jjp s; p = 0 

o 

[[G(r,c;r 0 )]] r=r o 
o 

where v is Poisson's ratio. And, of course, Green's function must still satis- 
fy the mixed boundary conditions (tb). 

Since r D is now the design variable, equation (12) is replaced by 

6G G r (r,c;r 0 )6r 0 (21) 

o 

where 6 r Q denotes the infinitesimal shift in the loeation of vanishing stiff- 
ness. 

It is desired to obtain an explicit representat ion for oG, analogous to 
equation ( 17 ). Toward this end, it is noted that 6G satisfies L 6 G =0 and there- 
fore 

(S(r6G rr + v6G r )) rr - (SCr^fiGp + v6G rr )) r = 0 ( 22 ) 

for a<r<r 0 and r 0 <r<b. Also 6G satisfies the same boundary conditions at a and 
b as does G. 

Before considering the jump conditions for 6G at r»r Cf some notational 
simplication can be obtained by noting that G(r,£;r 0 ) represents the response 
at the circle of radius r due to a unit load distributed along the circle of 
radius 5. Accordingly, let M (r,^;r 0 ) and Q (r,c;r 0 ) denote, respectively, the 
radial bending moment and shear force per unit length along the radius r due to 
the same unit load acting along the radius 5. Thus equations ( 20 ) simplify to 


rM ( p , ^ ; Pq ) | p~ p — 0 

0 

(23a) 

t[Q(p,C;r 0 )]] r=r C 

0 

(23b) 

[[G(r,^ ;p 0 ) JJp^p^ - 0 

(23c) 

sake of completeness, it is noted that Fi 

and Q are related to G through 

rM = - S(pG + vG ) 

rr r 


Pw “ - p(S(pG + vG ) ) 

rr r 

(2H) 


+ S(G + vrG ) 
r rr 


For the varied design whose stiffness vanishes at r Q + 6r Q , the jump con 
ditions analogous to (23a), (23b) and (23c) are, respectively, 
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rM (r,c;r 0 ) + r6K(r,6;r 0 ) | r=r ffir = 0 

o o 

(25a) 

[[rQ(r ,c;r 0 ) + r6Q(r,c;r 0 )]] r=r +6r = 0 

(25b) 

0 o 

C[G(r,c;r 0 ) + 6G(r,c;r 0 ) J] r=r +6p = o 

(25c) 


o o 


By expanding rM| r= ^ +6 ^, rQ| r=r + ^ and G| r=r + 6r in Taylor series about 
r^r 0 , and simplifying the result usfng equations ft5af, (25b) and (25c), the 
jump conditions 


r6K| r= p ^ - (rM)p| r=r 6r 0 

(26a) 

O ' 0 

~ 0 

(26b) 

0 o 

[ [ 6G] ] r=r - * CLGpillp^p 6r c 

(26c) 


are readily obtained. The quantities 6M and 6Q are implicity defined through 
equations (24). Thus 

r5M = - S(r6G rr + v6G r ) 

r6Q = - r(S(roG rr + v6G r )) r (27) 


+ S(6G r + vr6G rr ) 

The sensitivity 6G may now be determined explicitly by multiplying both 
sides of equation (22) by G and integrating the result from r-a to r=b. Thus 

6GU.O = + [[rQ(r,c)6G(r,0]] r = r 

o 

- [ [ ri-i ( r , q ) 6 G r ( r , O ] J r- r Q 

+ [[G r (r,c) rcfl(r,6) ]] r ^ r 

o 

- [[G(r,c)r 6 Q(r,C)]] r=r (28) 

o 

The second and fourth terms on the right hand-side of equation (28) vanish ty 
virtue of equation (23s) and equations (23c) and (26b), respectively. Moreover, 
equations (26c) and (26b) transform the first term into 


— rQ(r ,c)| r=r ^t-Gp(r,6)33p=p STq 

o o 

while the third term, obtained from equation (26a) and the equilibrium 
equation, becomes 


- rQ(r,0 | r=r “[G r (r,c)]] r= , 


lp=p LLUpU ,l,;jjp=p 

O O 


6r, 


Therefore, the design derivative 3C/3r 0 is given by 
r o)/9r 0 . ir 0 Q(p 0>t)[[Cp(rf6)]] 


r*r 


+ r 0 Q(r 0 ,C)[[G r (r ,^)]] r=r } 

o 


(29) 
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APPLICATIONS TO OPTIMAL DESIGN 


The usual method of obtaining structural optimality criteria associated 
with specific cost functionals relies on developing an appropriate variational 
formulation of the field equations (la) and (lb). Moreover, each cost 
functional requires a different variational formulation. In contrast to the 
historical approaches, the design derivatives specified by equations (10) and 
(17) or (29) can be used directly to determine the optimality criteria associ- 
ated with any cost functional without the need of a variational formulation. 

In order to illustrate the foregoing claim, structural optimality criteria will 
be derived for two different cost functionals: minimum response and minimum 

compliance. 

The optimality criteria associated with the design of a fixed-weight 
structure for minimum response at a specified location is considered first. 
Boundary conditions and loads are assumed known. It is desired to obtain a 
complete description of the design variable S including its singular points 
(locations of vanishing stiffness). Shield and Prager (ref. 5) obtained the 
optimality condition for this problem only after discovering the principle of 
stationary mutual potential energy. They did not address the question of locat- 
ing the singular points. However, at least one author (ref. 6) incorrectly 
assumed that such points can be obtained by requiring the response to be con- 
tinuously differentiable everywhere. 

Attention is now directed toward equation (la) and (lb) with g » h = 0. 
According to equation (3), the solution for the response is 

“ ■ 

Let the location y be specified and u(y) be a minimum. Thus 

u(y ) = ( f ( • ) , G( • ,y ) )q (30) 

whence 

6u(y) - (f( •), 6G(.,y)) fi (31) 

For the moment, it will be assumed that S is not singular anywhere. After sub- 
stituting equation (10) into equation (31) and changing the order of the resul- 
ting double integration, equation (31) becomes 

3E 

fiu(y) - - (Tu, 6STG(.,y)^ (32) 

The volume constraint may be easily handled through a Lagrange multiplier. Let 
v(S) denote the specific volume and A a Lagrange multiplier. Then the con- 
dition 6u(y) = C for all designs consistent with the constant volume constraint 
requires that the augmented functional 

- (Tu, ff 6STG <-,y)) fl + A (i.fj 6S)a = 0 

for all variations oS. Thus the optimality condition 

Tu.|| • TG(.,y) = A|| (33) 
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follows immediately. For beams, the equivalent to equation (33) was obtained 
in reference 5. The constant A appearing in equation (33) can be determined 
from the fixed-volume constraint. 

It was stated earlier that, in many applications, there can be no solution 
to the optimality and field equations unless a singular point occurs within the 
structure. In this case it is not possible to determine the location of the 
singularity from the optimality condition and field equations. This location 
must be considered on additional design variable, and consequently its location 
will be determined from an additional optimality condition. 

For simplicity, it will be assumed that the structure is a beam. In this 
case, equation (17) is substituted into equation (31). and the resulting double 
integrals are evaluated by reversing the order of the integrations. Thus 

6u(y) — — Q(x 0 )[[G z (z,y) ]] z - x ^Xq 

o 

- Q(x 0 ,y)[[u z (z)]] z=x 6x 0 (3*0 

o 

Since the specific volume v is independent of x D , the optimality condition to 
determine x 0 is obtained directly from equation (3 1 *). Thus 

Q(x 0 ) [ [G z (z ,y ) ]3 x=Xq 

+ Q(x 0 y)[[u z (z)]] x=Xo = 0 (35) 

Next, consider the problem of minimizing the compliance of a structure. 

The compliance C is defined to be the work done by the external loads. Thus 

C - (u,f) fl 


whence 


6C =■= (6u,f) n (36) 

Substitution of equation (32) into equation (36) yields 
/ 9E 

6C = - (Tu, gj 6STu) a 

Consequently, the optimality conditon for prescribed volume becomes 

Tu. ||- • Tu = A (37) 

Equation (37), in its various specific forms, has been derived by many authors 
for specific structures. In virtually all instances, the principle of minimum 
potential energy has been an ingredient necessary to the derivation. 

The location of any singular points may be determined in the same way that 
it was done for the minimum response design. In this case, equation (3^) is 
substituted into equation (36) to yield 

6C = - 2Q(x 0 )[[u z (z)]] z=x 6x 0 
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( 38 ) 


The optimality condition to determine x D , therefore, is 

[[u z (z)]] z ~ x = 0 

o 

As a final remark, it is pointed out, without elaboration, that the 
approach taken in this paper is easily generalized to transient structures. 
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